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ABSTRACT 

Context. The presence of non-circular and off-plane gas motions is frequently observed in the inner regions of disk galaxies. 
Aims. We have measured with integral-field spectroscopy the surface-brightness distribution and kinematics of the ionized gas in 
NGC 2855 and NGC 7049. These two early-type spiral galaxies were selected to possibly host a kinematically-decoupled gaseous 
component in orthogonal rotation with respect to the galaxy disk. 

Methods. We have modeled the ionized-gas kinematics and distribution of both galaxies assuming that the gaseous component is 
distributed either on two orthogonally-rotating disks or in a single and strongly warped disk. 

Results. In both galaxies the velocity field and distribution of the inner gas are consistent with the presence of an inner polar disk. In 
NGC 2855 it correponds to the innermost and strongly warped portion of the main disk. In NGC 7049 it is a central and geometrically- 
decoupled disk, which is nested in the main disk. 
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1. Introduction 

In a recent paper we pointed that about 50% of bright unbarred 
galaxies show a remarkable gas velocity gradient along the op- 
tical minor axis (Coccato et al. 2004). This phenomenon is ob- 
served all along the sequence of disk galaxies from SO's to Sm's. 
Since minor-axis velocity gradients are unexpected if the gas is 
moving onto circular orbits in a disk coplanar to the stellar one, 
we concluded that non-circular and off-plane gas motions are 
not rare in the inner regions of disk galaxies. Further support to 
this picture is given by the analysis of the velocity fields of large 
samples of lenticular galaxies (Sarzi et al. 2006), early (Falcon- 
Barroso et al. 2006) and late-type spirals (Ganda et al. 2006) 
observed with integral-field spectroscopy. 

The presence of a velocity gradient along the minor axis is 
characteristic of galaxies hosting an inner polar disk (IPD here- 
after). IPDs are small disks of gas and/or stars (R * 300 pc), 
which are located in the center of lenticular and spiral galax- 
ies and are rotating in a plane perpendicular to that of the main 
disk of their host. We call them IPDs instead of inner polar 
rings (Sil'chenko 2006) since there is no clearcut evidence they 
have an anular structure. Most of these orthogonally-rotating 
disks have been discovered in last few years (Corsini et al. 
2003; Sil'chenko & Afanasiev 2004; Shalyapina et al. 2004; 
Sil'chenko & Moiseev 2006; Coccato et al. 2005). The acquisi- 
tion of external gas via merging or accretion on nearly polar or- 
bits by a pre-existing galaxy, and the transfer of gas onto highly- 
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inclined anomalous orbits of a triaxial bulge or a bar, which is 
tumbling about its short axis, are both viable mechanism to build 
a orthogonally-rotating disk. According to these different sce- 
narios, it can be either a geometrically-decoupled structure or 
the inner portion of a strongly warped and larger gaseous disk, 
respectively. Therefore, constraining the structural properties of 
a sample of IPDs will give the clues to understand the processes 
driving their formation (see Sil'chenko 2006 for a review). 

Over the course of the last few years, we have undertaken 
a program aimed at detecting IPDs in disk galaxies using long- 
slit spectroscopic observations (Bertola et al. 1999; Sarzi et al. 
2000; Corsini et al. 2002, 2003; Coccato et al. 2004, 2005). In 
long-slit spectra the kinematic signature of an IPD is the pres- 
ence of a central velocity gradient along the disk minor axis and 
a zero-velocity plateau along the disk major axis, respectively. 
On the contrary, the presence of a gas velocity gradient along 
both major and minor axis is characteristic of gas moving onto 
elliptical orbits in the disk plane of a triaxial bulge (de Zeeuw & 
Franx 1989; Corsini et al. 2003). Our goal was to identify galax- 
ies, which are good candidates to host a IPD, to be followed up 
with integral-field spectroscopy at high spatial resolution in or- 
der to determine the size and orientation of such a kinematically- 
decoupled component. 

In this paper we present the analysis of the bidimensional 
kinematics and distribution of the ionized-gas component of the 
two most promising objects we found in our investigation. They 
are the nearby early-type spirals NGC 2855 and NGC 7049. For 
both galaxies the presence of a major-axis velocity plateau to- 
gether with a minor-axis velocity gradient is suggestive of an 
IPD as discussed by Corsini et al. (2002, 2003). The reader is 
refereed to them for a summary of the main properties of the 
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Table 1. Instrumental setup 
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Parameters 




Channels 






1 


2 


3 


4 


Grism 


HR_red 


HRjred 


HRj-ed 


HR_orange 


Order sorting filter 


GG475 


GG475 


GG475 


GG435 


Spatial resolution (" fiber -1 ) 


0.67 


0.67 


0.67 


0.67 


Reciprocal dispersion (A pixel -1 ) 


0.58 


0.58 


0.58 


0.62 


Readout noise (e - ) 


4.5 


3.2 


3.2 


3.7 


Gain(e - ADIT 1 ) 


1.9 


1.8 


1.7 


2.0 


Wavelength range (A) 


6300 - 8700 


6300 - 8700 


6300 - 8700 


5200 - 7600 


Instrumental FWHM at H« (km s -1 ) 


68 


68 


68 


70 



Table 2. Observing log 



Galaxy 


Date 


Observing block 


Exposure time 


NGC 2855 


27 Feb 2004 


154289 


2 x 19 min 


NGC 2855 


27 Feb 2004 


154292 


2x 19 min 


NGC 7049 


13 May 2004 


154293 


2x19 min 


NGC 7049 


14 May 2004 


154296 


2x19 min 



Notes - During the execution of the observing block 154293 (13 
May 2004) the data of channel 4 were lost because the Grism 
Exchange Unit of channel 4 was out of order. 



two galaxies. The paper is organized as follows. The integral- 
field spectroscopic observations and data reduction are described 
in Sect. [2] We present and model the ionized-gas kinematics and 
distribution in Sect. [3] and [4] respectively. Our conclusions are 
discussed in Sect. [5] 

2. Observations and data reduction 

The integral-field spectroscopic observations of NGC 2855 and 
NGC 7049 were carried out with the Very Large Telescope 
(VLT) at the European Southern Observatory (ESO) in Paranal 
(Chile) on February 26-27, 2004 and May 13-14, 2004. The 
Unit Telescope 3 (Melipal) mounted the Visible Multi Object 
Spectrograph (VIMOS) in the Integral Field Unit (IFU) con- 
figuration. The field of view of the four VIMOS channels was 
27" x 27" and it was projected onto a microlenses array. This 
was coupled to optical fibers which were rearranged on a linear 
set of microlenses to produce an entrance pseudoslit to the spec- 
trograph. The pseudoslit was 0795 wide and generated a total 
of 1600 spectra covering the field of view with a spatial reso- 
lution of 0767 per fiber. Each channel was equipped with either 
the HR_red or HR_orange high resolution grism and a thinned 
and back illuminated EEV44 CCD with 2048 x 4096 pixels of 
15x15 //m 2 . The details of the instrumental setup are given in 
Table Q] 

For each galaxy we obtained 4 x 19-min exposures. They 
were taken in service mode by executing two different observ- 
ing blocks of two exposures each, as listed in Table [2] Different 
quarz and arc lamp spectra were taken after every object expo- 
sure to ensure accurate flatfield correction and wavelength cali- 
bration, respectively. During observations seeing FWHM ranged 
between 078 and 170 as measured by the ESO Differential Image 
Meteo Monitor. 

For each VIMOS channel all the spectra were traced, iden- 
tified, bias subracted, flatfield corrected, corrected for relative 
fiber transmission, and wavelength calibrated using the routines 



of the ESO Recipe Execution pipelin^J Cosmic rays and bad 
pixels were identified and cleaned using standard MIDAS 1 rou- 
tines. We checked that the wavelength rebinning was done prop- 
erly by measuring the difference between the measured and pre- 
dicted wavelength for the brightest night-sky emission lines in 
the observed spectral ranges (Osterbrock et al. 1996). The re- 
sulting accuracy in the wavelength calibration is better than 5 
km s _1 . The intensity of the night-sky emission lines was used 
to correct for the different relative transmission of the VIMOS 
channels. The contribution of the night-sky emission lines was 
determined from a number of spectra by fitting a Gaussian to 
each line and a straight line to the adjacent continuum. We se- 
lected spectra where the night-sky emission lines did not over- 
lap with the [N II] /L16548, 6583 and Ha emission lines of the 
galaxy. Then night-sky emission lines were subtracted from all 
the available spectra. It was not possible to subtract the contribu- 
tion of the night-sky continuum alone since both galaxies cov- 
ered the entire field of view of the integral field unit. The pro- 
cessed spectra were organized in a datacube using the tabulated 
correspondence between each fiber and its position in the field 
of view. From the two exposures of each observing block we 
built a single datacube. In fact, the two spectra taken by the same 
fiber were coadded after checking that there was no pointing off- 
set between the two exposures. This was done by comparing the 
position of the intensity peaks of the two reconstructed images 
obtained by collapsing the datacubes along the wavelength di- 
rection. Finally, for each galaxy we coadded the two available 
datacubes by using the intensity peaks of the flux maps as a ref- 
erence for the alignment. In this way we produced a single dat- 
acube to be analyzed in order to derive the surface brightness 
and kinematic maps. 

3. Kinematics and distribution of the ionized gas 

3.1. Analysis 

In the coadded spectrum of the galaxy datacube we measured 
the position, full width at half maximum (FWHM), and uncal- 
ibrated flux of the Ha and [N II] -16583 emission lines by fit- 
ting a Gaussian to each line and a straight line to the adjacent 
continuum. This includes the contributions of both galaxy and 
night sky. No double-peaked emission line was observed. We 
neglected the [N II] /16548 line because it was not clearly de- 
tected in all the spectra, and we averaged the spectra of adja- 
cent fibers in the outer regions where the intensity of the Ho- and 
[N II] /16583 lines was low. We performed a 2 x 2 binning in 
order to reach a minimum signal-to-noise ratio S/N - 10. To 

1 ESOREX and MIDAS are developed and maintained by the European 
Southern Observatory 
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g. 1. Maps of the ionized-gas kinematics of NGC 2855 deriv 
)m the Ha (top panels) and [N II] /16583 (bottom panels) em 
sion lines. The field of view of each panel is 27" X 27". In each 
panel the spatial coordinates with respect to the galaxy center 
are given in arcsec. East is up and North is right. The ranges 
are indicated at the bottom of each panel. The white color is 
used for bins where no measurement was obtained on account of 
a signal-to-noise ratio below the selected thresold due to either 
faint emission line or poor fiber transmission. Left panel: Surface 
brightness in arbitrary linear units. Central panel: Heliocentric 
line-of-sight velocity in km s _1 without applying any correction 
for galaxy inclination. Isovelocity contours correspond to ve- 
locities after the subtraction of the systemic velocity. The zero- 
velocity contour is highlighted with a thicker solid line. Right 
panel: Line-of-sight velocity dispersion in km s _1 corrected for 
instrumental velocity dispersion. 



assess the data quality we show some examples of the profile of 
the [N II] /16583 emission line and its Gaussian fit in different 
regions of the two galaxies. They are plotted in Figs.[3]and|4]for 
NGC 2855 and NGC 7049, respectively. The central wavelength 
of the fitting Gaussian was converted into line-of-sight velocity 
in the optical convention. The standard heliocentric correction 
was applied. The Gaussian FWHM was corrected for the instru- 
mental FWHM and then converted into the intrisic line-of-sight 
velocity dispersion. The flux of the fitting Gaussian was con- 
verted into surface brightness according the fiber area. No flux 
calibration was performed. 

The result is shown in Figs. Q] and [2] for NGC 2855 and 
NGC 7049, respectively. 

To study the distribution of the ionized gas we fitted elliptical 
isophotes to the [N II] -16583 surface brightness map by means 
of the IRAI0 taks ELLIPSE. We first fitted ellipses allowing the 
center to vary to test the regularity of the gas distribution. We 
found no evidence of a varying center with the fits for NGC 
2885. The ellipse fits were then repeated with the ellipse cen- 
ter fixed. This was not the case of NGC 7049. The position of 
center was constant within the errors only for r < 4", while it 
changed by w 2" at larger radii. The radial profiles of the sur- 
face brightness, ellipticity and position angle derived from the 



2 IRAF is distributed by National Optical Astronomy Observatories, 
which is operated by Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National Science 
Foundation. 
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Fig. 2. Same as in Fig. Q]but for NGC 7049. 

NGC 2655 

I i i i I i i i I i i i I i i i I 




3 4 6 B 10 

r [arcsec] 



Fig. 5. Radial profiles obtained from the isophotal analysis of the 
[N II] /16583 surface-brightness map measured for NGC 2855. 
Top panel: Surface brightness calculated adopting an arbitrary 
zero point of 20 mag arcsec -1 . Middle panel: Ellipticity. Bottom 
panel: Position angle. In each panel the solid line corresponds to 
the prediction of the model which assumes that the ionized-gas 
component is distributed onto two orthogonal disks (see Sect. 
14. 2\ , while the dashed line corresponds to the prediction for the 
single disk model (see Sect. 14. 11 1. 



[N II] -16583 surface-brightness map are shown in Figs. [5] and [6] 
for NGC 2855 and NGC 7049, respectively. 

3.2. Results 

In NGC 2855 the ionized gas is concentrated in the nucleus and 
its distribution at larger radii follows the spiral pattern of the 
galaxy (Fig. [TJ. The abrupt change of ellipticity (from 0.05 to 
0.45) and position angle (from 105° to 155°) of the ellipses fit- 
ting the [N II] -16583 surface-brightness map (Fig. [5]) are the pho- 
tometric signatures of such a distribution of the gaseous compo- 
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Fig. 3. Central panel: Maps of the surface brightness and line-of-sight velocity field derived from the [N II] /16583 emission line for 
NGC 2855. The field of view, orientation, ranges and isovelocity contours are as in Fig.Q] Outer panels: The [N II] /16583 emission 
line profile (solid line) and its Gaussian fit (dashed line) in different positions of the field of view. Line profiles are normalized to 
the peak intensity to allow comparison. 
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Fig. 6. Same as in Fig.|5]but for NGC 7049. 

nent. The gas velocity field is complex. It is characterized by a 
S-shaped zero velocity line in the nuclear region. In the outer 



region it shows a regular rotation around the galaxy minor axis. 
The velocity dispersion increases from the outer parts towards 
the center. It peaks to »120 km s in Ha and to « 150 km s 
in [N II] /16583 and is systematically smaller for Ha with respect 
to [Nil] ^6583 (Fig. 

In NGC 7049 the gas is mostly concentrated in the nucleus 
and in a ring-like structure (Fig. This corresponds to the 
prominent dust lane, which is visible in the optical images of the 
galaxy (Sandage & Bedke 1994, Panel 74). The photometric pro- 
files derived from the [N II] /i6583 surface-brightness map (Fig. 
ID allows to constrain the different shape and orientation of the 
nuclear (e ~ 0.1, PA* 85°) and ring (e « 0.3, PA« 55°) regions. 
The gas velocity field of NGC 7049 is similar to that of NGC 
2855, displaying a S-shaped distortion in the nucleus and regular 
rotation around the galaxy minor axis at large radii. The velocity 
dispersion is systematically lower in Ha than in [N II] /16583. In 
the nuclear region shows a central dip surrounded by a double 
peak where it reaches w 120 km s~' in Ha and ss 170 km s~' in 
[N II] /I6583. Outwards it decreases smoothly with radius (Fig. 

EJ. 

In both galaxies the line-of-sight velocities measured from 
Ha and [N II] /16583 are consistent within the measurement er- 
rors. For this reason we modeled only the velocity field and 
surface-brightness distribution derived from the [N II] /i6583 
line. Typical errors on the line-of-sight velocity and velocity dis- 
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Fig. 4. Same as Fig.fjbut for NGC 7049. 



persions were m 10 and 20 km s l , including the uncertainties 
due to wavelength calibration. 

4. Analysis of the velocity field of the ionized gas 

The field of view of VIMOS/IFU allows us to map the velocity 
field of the kinematically-decoupled gaseous component, as well 
as that of the gas residing in the main disk of the galaxy. As a 
first step, we fit the whole velocity field by assuming that the gas 
is moving onto circular and coplanar orbits (Sec. 14.11 ). To this 
aim we masked the central region, where the S-shaped distorsion 
of the isovelocity countours is observed. We adopted the results 
of the circular velocity model to constrain the orientation and 
inclination of the main disk of the galaxy. Then we modeled the 
complex gas kinematics and surface-brightness distribution of 
the inner regions. We assumed that the gaseous component was 
distributed either on two distinct and orthogonally-rotating disks 
(Sec. 14.21 ) or onto a single and strongly warped disk (Sec. 14. 3b . 

In our analysis we fitted iteratively a model to the ob- 
served velocity field and surface-brightness distribution of the 
[N II] /16583 line using a non-linear least-squares minimiza- 
tion method. It was based on the robust Levenberg-Marquardt 
method implemented by More et al. (1980). The actual compu- 
tation was done using the MPFIT algorithm implemented by C. 
B. Markwardt under the IDL environment During each itera- 
tion the model was calculated on a subsampled pixel grid with 

3 The updated version of this code is available on 
http://cow.physics.wisc.edu~craigm/idl/idl.html 



the bin size of 07335 x 07335 (i.e., a subsampling factor of 2 x 2 
relative to the VIMOS pixel scale). This resulted the best com- 
promise between good sampling and reasonable computational 
time. The model was rebinned on a spatial scale of 0767 x 0767 
and the seeing effects were taken into account by convolv- 
ing the model with a Gaussian kernel with FWHM = 1". 
Finally, the model was compared to the observed velocity field 
and surface-brightness distribution. 

4.1. Single rotating disk 
4.1.1. Model calculation 

The model of the gas velocity field is generated assuming that 
the ionized-gas component is moving onto circular orbits in an 
infinitesimally thin disk with a negligible velocity dispersion. 
Hereafter this will be considered as the main disk of the galaxy. 

Let (X, Y, Z) be Cartesian coordinates with the origin in the 
center of the gaseous disk, the T-axis aligned along the line of 
nodes, and the disk plane confined to the {X, Y) plane. We can 
consider the gas motion in this plane by introducing the polar 
coordinates (R, (p) with origin in the center of the gaseous disk. 
It is R = y/X 2 + Y 2 (or R = y/(X - X c ) 2 + (Y - Y c ) 2 if the coor- 
dinates of the kinematic and geometric center of the gaseous disk 
are different) and cos <p = Y/R with <p counted counter-clockwise 
from T-axis. We assume that the gas circular velocity V at a 
given radius R is 

2 R 

V(R)= -V^arctan— (1) 
n R h 



6 



Coccato et al.: Inner polar disks in NGC 2855 and NGC 7049 



where V max and R/, are the maximal velocity and a scale radius, 
respectively. This empirical function has been proved to repro- 
duce fairly well the shape of optical rotation curves with the 
smallest number of free parameters (Courteau 1997). 

We now project the velocity field of the gaseous disk on the 
plane of the sky. Let (x, y, z) be Cartesian coordinates with the 
origin in the center of the gaseous disk, the x-axis pointing east- 
wards, the y-axis pointing northwards, and the z-axis along the 
line of sight directed toward the observer. The sky plane is con- 
fined to the (x, y) plane. In the reference frame of the sky (x, y, z), 
the gaseous disk is inclined by the zenithal and azimuthal angle i 
and 9, respectively. These angles correspond to the inclination of 
the disk (with i - 0° corresponding to the face-on case) and posi- 
tion angle of its line of nodes (with 9 counted counter-clockwise 
from the y-axis and 9-0° corresponding to the case with the 
line of nodes pointing northward). The transformation of the co- 
ordinates of the main disk from its reference frame (Z = 0) to 
the reference frame of the sky (x, y, z) is given by 



(2) 



' X" 




( X ) 


y 


= 91 


Y 


,z ) 




lo J 



where 



91 = 



cos 9 cos i sin 9 cos 9 sin i \ 
- sin 9 cos i cos 9 - sin 9 sin i 
- sin i cos i 



which leads to the usual equations in the sky plane (x,y): 

x — X cos 9 cos i + Y sin i 
y = -X sin cos ; + Y cos 9 

The ionized-gas velocity measured along the line of sight at 
a given sky point (x, y) is 



v(x, y) = V (R) sin i cos <p + V s , 



(4) 
is the 



where the circular velocity V is given in Eq. Q] and V, 
systemic velocity of the galaxy. 

The parameters of our model are the maximal velocity Vmax, 
scale radius Rh, inclination i and position angle 9 of the gaseous 
disk, and the systemic velocity V sys of the galaxy. The best-fitting 
model was obtained by masking the region of the observed ve- 
locity field characterized by the S-shaped distorsion of the isove- 
locity countours. The area of the masked region was 9x11 pixels 
and 7 x 8 pixels in NGC 2855 and NGC 7049, respectively. 

4.1.2. Results 

The velocity fields of the best-fitting disk models are compared 
to observed ones in Fig. [7] Best-fitting parameters are given in 
Tabled 

For the gaseous disk of NGC 2855, the position angle and 
inclination measured at large radii with the isophotal analysis 
(Sec. 13 - 1 b are not consistent with that derived from the kine- 
matic model (Tab. [3]). We explain such a discrepancy as due to 
the presence of the patchy spiral arm pattern, which character- 
izes the surface brightness distribution of the [N II] /16583 line 
at r > 4". This makes the photometric parameters not reliable 
to describe the geometry of the disk. The coordinates of the disk 
centers derived in the photometric and kinematic analysis (Tab. 
[3]) coincide within a pixel. 

The position angle, inclination and center derived from the 
photometric and kinematic analysis of the gaseous disk of NGC 
7049 are consistent within the errors. 
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Fig. 7. Model of the velocity field of NGC 2855 (upper panels) 
and NGC 7049 (lower panels) with a single rotating disk. The 
field of view, orientation, ranges and isovelocity contours are 
as in Fig.[T] Left panel: Observed velocity field. Central panel: 
Model. Right panel: Residuals. 



(3) Table 3. Model with a single rotating disk 



Parameter 


NGC 2855 


NGC 7049 


X c [pixel] 


-0.4±0.5 


-0.2±0.3 


Y c [pixel] 


0.4±0.6 


-0.6±0.3 


V sys [km s -1 ] 


1893+6 


2230±5 


en 


116.6+0.8 


58.1±0.7 


«'["] 


42±4 


60±2 


Vmax [km S~'] 


-344±22 


-459± 13 


Rh ["] 


3.9+0.7 


3.5+0.4 



In the rest of the paper we adopt as position angle and incli- 
nation of the main gaseous disks of NGC 2855 and NGC 7049 
the values we derived from the kinematic model. 



4.2. Orthogonally-rotating disks 
4.2.1. Model calculation 

The model of the gas velocity field is generated assuming that 
there are two different ionized-gas components which are mov- 
ing onto circular orbits in two orthogonally-rotating disks. They 
are the main galactic disk studied in Sect. I4.ll and an inner or- 
thogonal disk. Both disks have the same center, are infinitesi- 
mally thin, and have a negligible velocity dispersion. 

Let (X',Y',Z') be Cartesian coordinates with the origin in 
the center of the orthogonal disk, the T'-axis aligned along the 
line of nodes (defined as the intersection between the plane of the 
orthogonal disk and that of the main disk), and the plane of the 
orthogonal disk confined to the (X', Y') plane. In the reference 
frame of the main disk (X, Y, Z), the orthogonal disk is inclined 
by the zenithal and azimuthal angle 5 and y, respectively. These 
angles are the inclination of the inner disk (with 5 = 90° corre- 
sponding to the case of an orthogonal disk) and position angle 
of its line of nodes (with y counted counter-clockwise from the 
Y-axis and y — 0° corresponding to the case with the line of 
nodes of the orthogonal disk pointing in the same direction as 
the line of nodes of the main disk). A sketch of the orthogonal 
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disk in the reference frames (X, Y, Z) and (X', Y',Z') is shown in 
Figure [8] 

We assumed that the circular velocity of the gaseous compo- 
nent in the main disk and in the orthogonal disk are given by Eq. 
Q]and 



2 R' 
Vod(R') = -V^arctan— , 

71 K, 



(5) 



respectively. V' max and R', are the maximal velocity and scale ra- 
dius of the orthogonal disk. It is R' = VX' 2 + Y' 2 and cos (p' = 
Y' /R' with <p' counted counter-clockwise from T'-axis. 

The flux radial profile of the main gaseous disk was assumed 
to be exponential 



Fmd(R) = F + F x e 



-R/Rf 



(6) 



For the flux radial profile of the orthogonal disk we assumed a 
truncated exponential function 



F ( K'\-i F 'o + F 'i e ~ R ' IR ' f for R'<R' Q 
Fod(R)-{ q ° for R'>R' n 



(7) 



where R' Q corresponds to the radial extension of the orthogonal 
gaseous component. 

We now project the velocity field of the two gaseous disks 
on the plane of the sky. The transformation of the coordinates of 
the inner polar disk from its reference frame (X', Y',Z') to the 
reference frame of the main disk (X, Y, Z) is given by 



(8) 



(X s 
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where 
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' cosy cos 8 siny cosy sin 8 \ 
— sin y cos 8 cos y - sin y sin 8 
-sin 8 cos 8 



(9) 



The transformation of the coordinates of the main disk from 
its reference frame (X, Y, Z) to the reference frame of the sky 
(x, y, z) is given by 



(10) 
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where 



cos v cos i sint 
-sin 9 cos; cos I 
- sin i 



cos 9 sin i \ 
- sin 9 sin i 
cos i ) 



(11) 



The transformation of the coordinates of the orthogonal disk 
from its reference frame (X', Y',Z') to the reference frame of the 
sky (x, y, z) is given by 



(12) 



where 31 = SR25R1. 

The inclination i' of the orthogonal disk with respect to the 
sky plane is given by the direction cosine between the z-axis of 
the line of sight and the Z'-axis of the orthogonal disk 



' x) 




(X' \ 


y 


= 3? 


Y' 


y Z , 




\Z' ) 



The line of nodes of the orthogonal disk, defined as the in- 
tersection between the disk plane and the sky plane, is 

Z' = % x x + % 2 y + % 3 z = (14) 

where 5R' denotes the inverse of S JJ. In the sky plane it is 

% lX + % 2 y = (15) 

and therefore the position angle 9' of the apparent major axis 
of the orthogonal disk with respect to the y-axis is given by 
tanfl' =-9^ 2 /9^ r Thisis 

, sin #(sin 5 cos y cos i + cos 5 sin i) + cos 9 sin 6 sin y 

tan# = (16) 

cos 9(sm 6 cos y cos i + cos 5 sin i) - sin 9 sin 6 sin y 

The ionized-gas velocity measured along the line of sight at 
a given sky point (x, y) is 



v(x,y) = 
where 



VMp(x,y)fMp(x,y) + voD(x,y)fo D (x,y) 
f(x,y) 



VMo(x,y) = Vmd(R) sin; cos (f> + V sys 
fjm(x,y) = F M d{R)I cos i 



(17) 



(18) 
(19) 



are the line-of-sight velocity and surface brightness of the main 
disk, 



voD(x,y) = Vod(R) sin i' cos c 
foD(x,y) = F 0D (R)/ cos i' 



(20) 
(21) 



COS I 



9? 



3,3 



cos 6 cos i - sin 8 cos y sin i. 



(13) 



are the line-of-sight velocity and surface brightness of the or- 
thogonal disk, and f(x,y) - fiwD(x,y) + foo(x,y) is the observed 
surface brightness. 

For the main disk we adopted the inclination i and position 
angle 9 derived in Sect. 14. II and we assumed the inner disk to be 
orthogonal with respect to the main one within ±1°. Therefore 
the parameters of our model are V max , Rh, Fq, F\, and Rf for the 
main disk, V' mcu , R' h , F' Q , F[,R' F , R' Q , i' (with 89° < \i' - i\ < 91°) 
and 9' for the orthogonal disk, and the systemic velocity V SVJ of 
the galaxy. 

We started fitting a model only to the observed [N II] /16583 
surface-brighteness distribution (Fit #1 hereafter). Then we fit- 
ted a model only to the observed velocity field adopting for the 
main disk the maximal velocity V max and scale radius R/, derived 
in Sect. 14. 11 and for the orthogonal disk the inclination i' and po- 
sition angle 9' derived in Fit #1 (Fit #2 hereafter). Finally, we fit- 
ted simultaneously both the velocity field and surface-brightness 
distribution (Fit #3 hereafter). 

4.2.2. Results 

The surface-brightness distribution of the best-fitting models 
built in Fit #1 are shown in Fig. [9] The radial profiles of the 
surface brightness, ellipticity, and position angle derived from 
the model are compared to those from the isophotal analysis in 
Fig. and Fig. |6] for NGC 2855 and NGC 7049, respectively. 
The velocity field of the best-fitting models built in Fit #2 are 
shown in Fig.QT| The surface-brightness distribution and veloc- 
ity field of the best-fitting models built in Fit #3 are compared to 
observations in Fig.[T2land[T3l 

The surface-brightness distribution and velocity field of 
NGC 2885 can not be interpreted as due to presence of two dif- 
ferent orthogonally-rotating disks of gas. The largest deviations 
in the residual maps are observed in the nuclear region. 
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A y 




Fig. 8. Panel A: Orientation of the plane (X, Y) of the main 
disk of the galaxy with respect to the sky plane (x,y). Panel B: 
Orientation of the plane (X', Y') of the orthogonal disk with re- 
spect to the plane (X, Y) of the main disk. 



For NGC 7079 the surface-brightness distribution of the 
model built in Fit #1 is consistent with observations suggesting 
the presence of an inner gaseous component which is orthogonal 
with respect to the main gaseous disk of the galaxy. However, as 
in the case of NGC 2855, the observations are not reproduced 
by the models built in Fit #2 and #3. We conclude that the circu- 
lar velocity curve of the orthogonal disk differs from Eq. We 
derive the velocity field of the orthogonal disk as 




5 10 15 S 10 15 -3 

Arb. Units Arb Units Arb. Unit 



Fig. 9. Model of the surface-brightness distribution of NGC 
2855 (upper panels) and NGC 7049 (lower panels) with two 
orthogonally-rotating disks (Fit #1). The field of view, orienta- 
tion, ranges and isovelocity contours are as in Fig.Q] Left panel: 
Observed surface-brightness distribution. Central panel: Model. 
Right panel: Residuals. 



NGC 7049 [Kill 




Fig. 10. The velocity field of the orthogonal-rotating disk of 
NGC 7049, as obtained from Eq. [22] The field of view of the 
panel is 9" x 9". East is up and North is right. The range is indi- 
cated at the bottom of the panel. Isovelocity contours correspond 
to velocities after the subtraction of the systemic velocity. 

Table 4. Model with two orthogonally-rotating disks (Fit #1) for 
NGC 7049 



v(x, y)f(x, y) - v M d(x, y)fMD(x, y) 
vod = (22) 

JOD 



and show it in Fig. [10] We tested different parametric prescrip- 
tions (Brandt 1960; Freeman 1970; Bertola et al. 1991; Persic et 
al. 1996) to fit the velocity field assuming the ionized gas is on 
circular orbits on a plane. But none of them gave better results 
than that obtained adopting Eq. 5. This suggests that proba- 
bly the gas kinematics plotted in Fig. [10] is not characterized 
by circular motions. The limited radial extension of the re- 
gion does not allow us to adopt more complicated functions 
as well as a larger number of free parameters to take into 
account for non-circular motions. 



Parameter 


F 


0.01±0.01 


F\ 


1.6±0.2 


Rf ["] 


15+2 




0.01±0.01 




22±1 


R' F ["] 


1.13±0.08 


7 [deg] 


23±4 


S [deg] 


91±1 


*o ["] 


3.6+0.1 
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Fig. 11. Model of the velocity field of NGC 2855 (upper panels) 
and NGC 7049 (lower panels) with two orthogonally -rotating 
disks (Fit #2). The field of view, orientation, ranges and isove- 
locity contours are as in Fig. [T] Left panel: Observed velocity 
field. Central panel: Model. Right panel: Residuals. 
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Fig. 12. Model of the surface -brightness distribution {upper pan- 
els) and velocity field lower panels) of NGC 2855 with two 
orthogonally-rotating disks (Fit #3). The field of view, orienta- 
tion, ranges and isovelocity contours are as in Fig.Q] Left panel: 
Observed data. Central panel: Model. Right panel: Residuals. 



4.3. Warped disk 
4.3.1. Model calculation 

The model of the gas velocity field is generated approximating 
the warped gas distribution with a series of concentric and circu- 
lar wires. 

Let (X' n ', Y' n ',Z' n ') be Cartesian coordinates with the origin in 
the center of the n-th wire, the i^'-axis aligned along the line 
of nodes defined as the intersection between the plane of the 
wire (X' n ', Y") and that of the main disk (X, Y). In the reference 
frame of the sky (x, y, z), the wire is inclined by the zenithal and 
azimuthal angle 5„ and y„, respectively. 

The transformation of the coordinates of the n-th wire from 
its reference frame (X'^, Y" ,Z' t [) to the reference frame of the 




NGC. 7049 [Ml") Fit p ^i: 1'^ [HI] residua 



Arb. Units 
NGC | 7049 [NIIl .residua. 



a? 



2000 2200 2400 2600 2000 2200 2400 2600 -100 -50 50 100 

V [km s~*] V [km s^] V [km sr^ 

Fig. 13. Same as in Fig.[l2]but for NGC 7049. 

sky (x, y, z) is given by Eq. [T2] where S and y are replaced by 6„ 
and y„, respectively. Similarly, the inclination i% and the position 
angle 0% of the apparent major axis of the n-th wire are given 

by 

cos = cos 6 n cos i - sin 5 n cos y„ sin i. (23) 
and 

,, sin #(sin S„ cos y„ cos i + cos 6 n sin i) + cos 8 sin 6„ sin y„ 

tan#„ = 

cos 0(sin 6„ cos y„ cos i + cos 8„ sin z) - sin 6 sin 6„ sin y n 

We assumed that the circular velocity of the gaseous compo- 
nent in the warped disk and is given by 



(24) 



V WD «) = -V; , «arctan§ r> 
n h 



(25) 



where V" mK and R'J are the maximal velocity and scale radius of 

the warped disk. It is R'; = ^X"? + F," 2 and cos^,' = Y%/R% 
with <f>' n ' counted counter-clockwise from F ( "-axis. 

The ionized-gas velocity and flux measured along the line of 
sight at a given sky point (x, y) where M wires are observed are 
given respectively by 



v(x,y) 



and 



£„=i Vn(x,y)f„(x,y) 



f(x,y) = Mf n (x,y), 
where 

v„(x,y) = V WD (R'^) sin i" cos (/>',; + V„ 



(26) 



(27) 



(28) 



is the line-of-sight velocity and f„(x,y) is the line-of-sight flux 
of the n-th wire. Eq. [27] implies that all the wires observed at 
a given position on the sky give the same contribution to the 
observed surface brightness. This is supported by the fact that 
all the wires observed at a given position are supposed to have 
similar radii, warping and twisting. Moreover, this choice al- 
lowed to better reproduce the inhomogeneous distribution of 
surface brightness due to the spiral arms of NGC 2855 and ring- 
like structure of NGC 7049. The attempt to model the surface- 
brightness distribution adopting an exponential radial profile was 
not successful. 
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Table 5. Model with a warped disk for NGC 2885 



Parameter 



Vmax [km S '] 


-373±16 


R„ ["] 


4.3±0.6 


k, ["] 


3.4±0.1 


h ["] 


0.82±0.03 


h [deg] 


-79.0±2 


Co [deg] 


85.7±0.7 


d [deg/"] 


-10.11±0.07 


V IVS [km s '] 


1885±3 



We evaluated the warping and twisting of the gaseous com- 
ponent from the radial profiles of ellipticity and position angle 
obtained from the isophotal analysis of the surface-brightness 
map of the [N II] /16583 line in Sect. O We found that the fol- 
lowing empirical functions 



SniK) = 

and 

y n (R'£) = c + aR'n 



-k 2 j K ~ k o 
arctan 

7T \ h 



h 

2 



(29) 



(30) 



reproduced fairly well the radial profiles of 5„ and y n we derived 
from Eq. [23] and Eq. [24] by adopting for i and 6 the values of 
the main disk fitted in Sect. 14.11 This constraints the warped 
structure to be confined in the inner region of the galaxy. 

Therefore the parameters of our model are V" wx and R'/ for 
the circular velocity, ko, k\, and for warping, cq and c\ for 
twisting, and V sys for systemic velocity. Forty wires were used 
to build the warped structure and cover the whole velocity field 
with at least one wire for each pixel bin. 

4.3.2. Results 

The velocity fields of the best-fitting models and their residuals 
are are shown in Fig. [14] Best-fitting parameters are given in 
Table [5] 

The warped model is able to reproduce the S-shaped zero 
velocity line we observed in the nuclear region of NGC 2855 
with velocity residuals lower than a few tens of km s _1 . At large 
radii the position angle and inclination of the warped disk are 
consistent with those of the main disk of the galaxy. At small 
radii, the warped disk is almost orthogonal to the main disk (Fig. 

The velocity field of NGC 7049 can not be interpreted as 
due to presence of a single gaseous component distributed in a 
warped disk. In fact, the deviations of the residual map are larger 
than those we found for the models with both a single disk and 
two orthogonally-rotating disks. 



5. Conclusions 

We have measured with integral-field spectroscopy the surface- 
brightness distribution and kinematics of the ionized gas in 
NGC 2855 and NGC 7049. These two early-type spiral galaxies 
were selected to possibly host an IPD on the basis of the analysis 
of their long-slit spectra (Corsini et al. 2002, 2003). 

In NGC 2855 the ionized-gas distribution peaks in the nu- 
cleus and follows the patchy spiral pattern of the galaxy. In 
NGC 7049 the gas is mostly concentrated in the nucleus and 
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Fig. 14. Model of the velocity field of NGC 2855 (upper pan- 
els) and NGC 7049 (lower panels) with a warped disk. The field 
of view, orientation, ranges and isovelocity contours are as in 
Fig.[T] Left panel: Observed velocity field. Central panel: Model. 
Right panel: Residuals. 




T t [Nil] | residue 



Fig. 15. Model of the surface-brightness distribution of NGC 
2855 with a warped disk. The field of view, orientation, ranges 
and contours are as in Fig. Q] Left panel: Observed surface- 
brightness distribution. Central panel: Model. Right panel: 
Residuals. 



in a ring-like structure. The ionized-gas velocity field of both 
galaxies is characterized by a S-shaped line of zero velocity. It 
is aligned with the disk minor axis in the inner regions (r < 8") 
and with the disk major axis at larger radii. This remarkable mis- 
alignement between the kinematical and photometrical major 
axes gives rise to the central velocity gradient and zero-velocity 
plateau observed along the disk minor and major axis, respec- 
tively. These kinematical features can not be explained as due to 
non-circular gas motions in the principal plane of the disk caused 
by the triaxial potential of the bulge or a bar (de Zeeuw & Franx 
1989; Corsini et al. 2003). 

As a first step we modeled the velocity field, except for the 
inner region where the S-shaped distorsion of the isovelocity 
countours is observed. We assumed that the gas is moving onto 
circular orbits in an infinitesimally thin disk with a negligible 
velocity dispersion. This allowed us to constrain the orientation 
and inclination of the main disk of the galaxy. Then we modeled 
the ionized-gas kinematics and distribution of both the inner and 
outer regions. We assumed that the gaseous component is dis- 
tributed either on two orthogonally-rotating disks or in a sin- 
gle and strongly warped disk. In both galaxies the velocity field 
and distribution of the inner gas are consistent with the presence 
of an IPD, which is in orthogonal rotation with respect to the 
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outer gas. In NGC 2855 the IPD correponds to the innermost 
and strongly warped portion of the main disk. In NGC 7049 it is 
a central and geometrically-decoupled disk, which is nested in 
the main disk. 

We exclude that warped gaseous disk of NGC 2855 is due 
moving onto anomalous orbits in the triaxial bulge or in bar that 
is tumbling around the minor axis. In fact, we do not observe any 
bar structure neither in optical (Corsini et al. 2002) nor in near- 
infrared images (Peletier et al. 1999; Mollenhoff & Heidt 2001) 
of the galaxy. Moreover, there is no evidence of gas in retrograde 
motion relative to the stars at large radii from the galaxy center, 
where the gas orbits are expected to be highly inclined with re- 
spect to the figure rotation (van Albada et al. 1982; Friedli & 
Benz 1993). On the other hand, the galaxy is acquiring exter- 
nal material in a direction close to the disk minor axis (Malin & 
Hadley 1997). We suggest that the central portion of the warped 
disk is composed by accreted gas, which is settling in the plane 
perpendicular to the long axis of the bulge. It is in orthogonal ro- 
tation with respect to the outer gaseous component, which lies on 
the plane perpendicular to the long axis of the bulge. This is the 
one containing the main disk of the galaxy. In NGC 7049 the two 
orthogonally-decoupled gaseous components are already settled 
onto the principal planes of the triaxial bulge. 

To date the presence of an inner polar disk has been con- 
firmed by means of integral-field spectroscopy in more than 20 
early-type disk galaxies (Sil'chenko 2006 and this paper). In 
spite of the growing statistics, the question about their forma- 
tion is still open. Different qualitative arguments support either 
an external or internal origin. They are the acquisition of external 
gas via accretion on nearly polar orbits by the triaxial bulge of a 
pre-existing galaxy, and the transfer of gas onto highly-inclined 
anomalous orbits of a triaxial tumbling bulge or bar, respectively. 
The comparison between the observed kinematics and distribu- 
tion of gas in bulges and bars hosting an inner polar disk with 
gas dynamical models including external acquisition and inter- 
nal evolution in a variety of triaxial potentials are now required 
to address this issue. 
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